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ABSTRACT 


Self-Propelled  Semi-Submersibles  now  transport  an  estimated  75%  of  cocaine 
originating  from  Colombia  and  headed  for  the  United  States.  There  are  several  types  of 
search  platforms  (i.e.,  units  to  detect,  classify,  and  interdict)  being  employed  by  the  Joint 
Interagency  Task  Force  South  to  combat  the  semi-submersibles.  We  use  a  defender- 
attacker  optimization  model  to  maximize  the  defender’s  probability  of  successful 
detection  and  classification  of  the  semi-submersible  through  the  advantageous  disposition 
of  these  search  platforms  against  an  intelligent  attacker  operating  the  semi-submersible. 
We  assume  the  attacker  has  imperfect  knowledge  of  defender  platform  disposition  but  is 
aware  that  there  are  defenders  that  must  be  avoided.  Given  this  assumption,  the  solution 
to  the  defender-attacker  model  is  a  mixed  (i.e.,  probabilistic)  strategy  for  the  defender  and 
a  least-risk  path  for  the  attacker.  We  demonstrate  our  defender- attacker  model  with  both 
an  Eastern  Pacific  and  a  Caribbean  scenario  using  five  representative  search  platform 
types  whose  detection  and  classification  performance  vary  by  platform,  and  by 
geography.  In  each  of  these  cases,  we  find  that  our  model  prescribes  a  face-valid 
defensive  plan;  defenders  take  advantage  of  geography  by  positioning  at  chokepoints  in 
constrained  waterways,  and  they  provide  coverage  near  attacker  origins  and  destinations 
in  the  less  geographically-constrained  scenarios. 
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EXECUTIVE  SUMMARY 


The  Self  Propelled  Semi-Submersible  (SPSS)  is  not  a  true  submarine,  but  a 
surface  vessel  with  very  low  freeboard.  Maritime  drug  traffickers  in  the  Caribbean  and 
Eastern  Pacific  have  historically  used  a  combination  of  commercial  fishing  vessels  and 
smaller,  harder-to-detect  “go  fast”  boats.  The  SPSS  combines  the  low  profile  of  a  “go- 
fast”  boat  with  the  long  range  and  high  capacity  of  a  fishing  boat.  Given  these 
advantages  the  SPSS  has  become  an  increasingly  popular  choice  for  smugglers 
trafficking  drugs  from  Columbia  to  the  Pacific  Coast  of  Mexico.  SPSSs  now  account  for 
up  to  75%  of  the  estimated  600  tons  of  cocaine  originating  from  Colombia  annually.  It  is 
not  lost  on  us  that  an  SPSS  can  also  carry  passengers,  or  even  a  weapon  of  mass 
destruction. 

There  are  several  types  of  search  platforms  (i.e.,  units  to  detect,  classify,  and 
interdict  SPSSs)  being  employed  by  the  Joint  Interagency  Task  Force  (JIATF)  South  to 
combat  the  SPSS.  These  include  surface  ships  and  aircraft  from  several  partner  nations 
and  U.S.  agencies.  We  also  consider  employing  non-traditional  stationary  acoustic 
sensors.  Submarines  are  worth  considering  given  the  SPSSs  relatively  loud  acoustic 
signature.  Though  a  variety  of  platforms  that  can  be  used  against  the  SPSS,  there  are 
generally  only  a  few  available  at  any  given  time. 

Because  searching  for  a  single  SPSS  can  tie  up  many  search  assets,  current  SPSS 
interdiction  relies  heavily  on  actionable  prior  intelligence.  In  the  absence  of  prior 
intelligence  operational  commanders  rely  on  information  such  as  cones  of  courses  and 
speed  versus  time  with  consideration  for  background  effects  such  as  weather  to  help 
determine  the  most  likely  routes  for  the  SPSS.  There  is  no  computer-based  tool  to  aid 
these  operational  platform  employment  decisions. 

We  introduce  a  new  planning  aid  for  operational-level  mission  planning  of  SPSS 
search.  This  planning  aid  will  provide  optimal  placement  and  disposition  of  cooperating 
friendly  search  platforms  while  considering  the  intelligent  response  of  the  enemy  SPSS 
operators  to  search  efforts.  A  successful  operation  against  an  SPSS  requires  detection, 


xv 


classification,  and  interdiction.  The  planning  aid  addresses  detection  and  classification 
only,  which  are  arguably  the  two  most  difficult  stages  of  the  operation. 

We  use  a  defender- attacker  optimization  model  to  maximize  the  defender’s 
probability  of  successful  detection  and  classification  of  the  SPSS  through  the 
advantageous  disposition  of  these  search  platforms  against  an  intelligent  attacker 
operating  the  semi-submersible.  We  assume  the  attacker  has  imperfect  knowledge  of 
defender  platform  disposition  but  is  aware  that  there  are  defenders  that  must  be  avoided. 
Given  this  assumption,  the  solution  to  the  defender-attacker  model  is  a  mixed  (i.e., 
probabilistic)  strategy  for  the  defender  and  a  least-risk  path  for  the  attacker. 

We  demonstrate  our  defender-attacker  model  with  both  an  Eastern  Pacific  and  a 
Caribbean  scenario  using  five  representative  search  platform  types  whose  detection  and 
classification  performance  vary  by  platform,  and  by  geography.  In  each  of  these  cases, 
we  find  that  our  model  prescribes  a  face-valid  defensive  plan;  defenders  take  advantage 
of  geography  by  positioning  at  chokepoints  in  constrained  waterways,  and  they  provide 
coverage  near  attacker  origins  and  destinations  in  the  less  geographically-constrained 
scenarios. 


xvi 


ACKNOWLEDGMENTS 


First  and  foremost  I  have  to  thank  my  advisor  Distinguished  Professor  Jerry 
Brown  for  his  expertise  and  encouragement.  This  thesis  would  not  have  been  possible 
without  him.  I  would  like  to  thank  my  second  reader  Senior  Lecturer  Jeff  Kline  for 
helping  me  find  a  topic  that  both  used  my  professional  background  and  allowed  me  to 
pursue  my  interest  in  attacker-defender  models.  Many  thanks  go  to  Associate  Professor 
Matt  Carlyle  for  spending  countless  hours  working  with  me  on  the  code  when  the 
inevitable  bugs  reared  their  ugly  heads.  Finally  I  have  to  thank  my  whole  family  for 
allowing  me  to  devote  so  much  time  to  this  thesis.  My  incredible  wife  Elizabeth  took 
care  of  business  despite  working  on  her  own  very  demanding  curriculum  at  the  Defense 
Language  Institute.  My  son  Andrew  and  my  daughter  Cate  embraced  added 
responsibility  to  allow  their  parents  to  reach  their  academic  goals. 


XVII 


THIS  PAGE  INTENTIONALLY  LEFT  BLANK 


I.  INTRODUCTION 


This  thesis  develops  an  operational  level  mission-planning  aid  to  counter  Self 
Propelled  Semi-Submersibles  (SPSS)  currently  being  used  by  drug  smugglers  in 
Colombia  to  move  large  volumes  of  cocaine  toward  the  U.S.  This  planning  aid  will 
suggest  optimal  cooperative  placement  and  disposition  of  friendly  search  platforms  while 
considering  the  intelligent  response  of  the  enemy  SPSS  operators  to  our  search  efforts.  A 
successful  operation  against  an  SPSS  requires  detection,  classification,  and  interdiction. 
This  thesis  addresses  detection  and  classification  only,  which  are  arguably  the  two  most 
difficult  stages  of  the  operation. 

The  primary  use  of  this  research  is  detennining  resource  versus  risk  to  gauge  the 
benefit  of  additional  platforms  and  capabilities  or  the  effect  of  eliminating  platfonns  and 
capabilities.  The  goal  is  to  discover  the  best  mix  of  platforms  to  use  against  the  SPSS  as 
well  as  their  optimal  cooperative  disposition  in  terms  of  placement  and  the  spread  of  their 
search  effort.  Major  inputs  to  the  model  are  positions  where  the  SPSS  can  enter  the  sea 
space  and  the  goal  positions  were  it  can  leave  the  sea  space.  A  smaller  number  of  these 
positions  implies  better  friendly  intelligence  and  will  result  in  improved  search  platform 
performance.  By  increasing  or  decreasing  the  number  of  entry  and  goal  positions,  the 
value  of  intelligence  can  be  captured. 

A.  SPSS  THREAT 

The  SPSS  is  not  a  true  submarine,  but  a  surface  vessel  with  very  low  freeboard 

(Figures  1  and  2).  A  typical  SPSS  has  a  hull  that  only  rises  a  foot  above  the  waterline. 

Maritime  drug  traffickers  in  the  Caribbean  and  Eastern  Pacific  have  historically  used  a 

combination  of  commercial  fishing  vessels  and  smaller,  harder  to  detect  “go  fast”  boats. 

The  SPSS  combines  the  low  profile  of  a  “go-fast”  boat  with  the  long  range  and  high 

capacity  of  fishing  boat.  SPSSs  have  been  in  limited  use  for  nine  years,  but  they  have 

become  an  increasingly  popular  choice  for  Columbian  drug  smugglers  in  the  Eastern 

Pacific  because  efforts  at  policing  fishing  boats  have  improved  in  recent  years  (Bajak, 

2008).  Twenty-three  SPSSs  were  launched  from  Colombia  between  2006  and  2007 
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(Wilkenson,  2008).  In  2008,  nearly  70  were  seen  or  captured  (StrategyPage,  2009). 
SPSSs  now  convey  up  to  75%  of  the  estimated  600  tons  of  cocaine  originating  from 
Colombia  annually.  The  vessels  have  been  operating  nearly  exclusively  in  the  Pacific 
making  the  transit  from  Colombia  to  Mexico  or  Guatemala  (Brown,  2009).  SPSS 
technology  is  easily  exportable  to  other  locations.  At  least  one  has  been  used  by  a  local 
drug  gang  in  Spain  (StrategyPage  2009). 

Though  SPSSs  are  currently  being  used  exclusively  for  drug  trafficking,  they 
could  easily  be  used  to  transport  any  manner  of  contraband.  Both  the  Director  of  Joint 
Interagency  Task  Force  (JIATF)  South  and  Commander,  United  States  Southern 
Command  (COMUSSOUTHCOM)  have  suggested  that  SPSSs  could  be  used  to  support 
terrorism  (Brown,  2009).  COMUSSOUTHCOM  cited  the  delivery  of  weapons  of  mass 
destruction  as  a  worst  case  scenario  (USSOUTHCOM,  2008). 

SPSS  have  small  visual,  infra-red  and  radar  signatures  making  them  difficult  to 
detect  by  non-acoustic  sensors.  Recently  SPSS  designers  have  taken  particular  care  to 
reduce  infra-red  signature  through  the  use  of  exhaust  cooling  systems  (Brown,  2009). 
Classification  is  also  difficult  for  most  sensors.  Because  each  SPSS  is  purpose-built,  and 
configurations  change  constantly,  there  is  no  telltale  signature  for  an  SPSS.  They  look 
similar  to  sailboats  on  radar  and  sound  similar  to  fishing  boats  acoustically.  As  a  result, 
an  SPSS  can  easily  blend  into  legitimate  maritime  traffic.  These  difficulties  are 
compounded  by  a  vast  area  of  operation  twice  the  size  of  the  continental  U.S.  It  is  little 
wonder  that  only  an  estimated  10%  of  SPSS  are  interdicted  (Wilkenson,  2008).  Many  of 
these  successful  interdictions  are  the  result  of  intelligence  information  at  the  source 
(StrategyPage,  2009).  Our  odds  of  success  can  be  improved  either  through  improved 
intelligence  or  through  improved  maritime  domain  awareness.  The  focus  of  this  study  is 
on  improving  maritime  domain  awareness  through  the  optimized  employment  of  scarce 
search  assets,  but  the  value  of  intelligence  can  also  be  captured  as  a  byproduct  of  the 
analysis. 
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Figure  1.  Various  SPSS  Designs  at  Bahia  Malaga  Navy  Base  Colombia  (From  Bajak, 

2008). 


Figure  2.  Captured  SPSS  at  NAS  Key  West  (From  Brown,  2009). 
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B. 


COUNTERING  THE  SPSS  THREAT 


1.  Search  Platforms 

There  are  several  types  of  search  platforms  (i.e.,  units  to  detect,  classify,  and 
interdict  SPSSs)  being  employed  by  the  Joint  Interagency  Task  Force  (JIATF)  South  to 
combat  the  SPSS  (e.g.,  Figures  3-5).  These  include  ships  such  as  partner  nation  frigates, 
U.S.  Navy  frigates  and  U.S.  Coast  Guard  cutters.  Aircraft  can  include  helicopters 
embarked  from  these  ships  as  well  as  land-based  maritime  patrol  aircraft  (MPA)  and 
airborne  early  warning  (AEW)  aircraft  operated  by  partner  nations,  U.S.  Air  Force,  U.S. 
Navy,  U.S.  Coast  Guard,  and  U.S.  Customs  and  Border  Patrol.  Consideration  is  also 
being  given  to  employing  non-traditional  acoustic  sensors  such  as  Seaweb  (Figure  6),  a 
proposed  autonomous  network  of  acoustic  sensors  that  communicate  with  each  other 
through  acoustic  modems  (Rice,  2009).  Submarines  are  considered  given  the  SPSSs 
relatively  loud  acoustic  signature  (at  least  relative  to  what  submarines  traditionally  hunt). 
Though  a  variety  of  platforms  that  can  be  used  against  the  SPSS,  there  are  generally  only 
a  few  available  at  any  given  time. 


Figure  3.  U.S.  Coast  Guard  C-130  Maritime  Patrol  Aircraft  (From  Wikipedia.org  2009). 
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Figure  4.  U.S.  Navy  Frigate  (From  Wikipedia.org  2009). 


Figure  5.  U.S.  Customs  and  Border  Patrol  P-3  Airborne  Early  Warning  Aircraft  (From 

Wikipedia.org  2009). 
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Figure  6.  Seaweb  is  an  Autonomous  Networked  Array  of  Acoustic  Sensors  (From  Rice 

2008). 


2.  Current  Planning  Methods 

Because  searching  for  a  single  SPSS  can  tie  up  many  search  assets  current  SPSS 
interdiction  relies  heavily  on  actionable  prior  intelligence.  This  greatly  reduces  the 
search  area.  In  the  absence  of  prior  intelligence  operational  commanders  rely  on 
information  such  as  cones  of  courses  and  speed  versus  time  with  consideration  for 
background  effects  such  as  weather  to  help  determine  the  most  likely  routes  for  the  SPSS. 
Care  is  taken  to  not  be  predictable  in  employing  search  assets,  since  smugglers  have 
demonstrated  the  ability  to  avoid  known  drug  interdiction  efforts.  There  are  no 
computer-based  tools  to  aid  these  operational  platfonn  employment  decisions.  There  are 
tactical  level  computer  based  planning  aids  such  as  the  Personal  Computer-Based 
Interactive  Multisensor  Analysis  Training  System  (PCIMAT)  designed  for  anti¬ 
submarine  warfare  that  can  help  predict  the  performance  of  acoustic  sensors.  However, 
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these  decision  aids  only  predict  the  perfonnance  of  individual  sensors;  they  cannot 
coordinate  the  simultaneous,  cooperative  disposition  of  several  platfonns  in  an  optimal 
manner  against  an  intelligent  adversary. 
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II.  MODEL  DEVELOPMENT 


A.  NETWORK  REPRESENTATION 

We  model  the  maritime  environment  using  a  mesh  network  (Figure  7).  The  ocean 
area  is  broken  down  into  60nm  square  cells  with  a  node  in  the  middle  of  each.  Each  cell 
node  is  connected  by  an  arc  to  and  from  adjacent  cells  unless  there  is  an  obstruction  to 
navigation,  such  as  land.  SPSSs  typically  begin  their  journey  from  small  river  estuaries 
or  sparsely-populated  coastline  where  our  maritime  platforms  cannot  effectively  operate. 
Therefore,  the  entry  cells  are  only  navigable  by  the  SPSS  (hereafter  called  the 
“attacker”).  SPSSs  typically  deliver  their  cargo  offshore  where  maritime  search 
platforms  can  operate.  Therefore,  goal  cells  are  navigable  by  both  the  attacker  and 
defender. 
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Figure  7.  Example  Network  Representation.  Each  dot  represents  the  center  of  a  60nm 
square  cell  of  sea  space,  and  each  line  connecting  a  pair  of  dots  represents  an 
adjacency  indicating  it  is  possible  to  navigate  directly  between  these  cells. 


We  assume  the  attacker  behaves  intelligently,  and  maximizes  his  probability  of 
evading  detection  and  classification  by  choosing  a  directed  path  from  an  advantageous 
entry  cell  node  to  a  an  advantageous  goal  cell  node  that  has  the  maximum  joint 
probability  of  evading  detection  and  classification  along  all  arcs  in  the  path.  If  we 
assume  independence  between  arcs,  the  joint  probability  that  the  attacker  will  evade 
detection  and  classification  is  the  product  of  the  evasion  probabilities  over  each  of  the 
arcs  the  attacker  transits.  Maximizing  the  sum  of  the  logarithms  of  the  evasion 


9 


probabilities  is  equivalent  to  maximizing  the  product  of  these  probabilities  because  the 
logarithm  is  a  single-valued  increasing  function  (e.g.,  Abdul-Ghaffar,  2008). 

B.  MISSION  GENERATION 

The  first  step  in  determining  the  arc-based  log  evasion  probabilities  is  the 
generation  of  all  possible  defender  missions  and  their  pre-calculated  cell-based  log 
evasion  probabilities.  A  mission  has  two  parameters,  the  set  of  cells  searched, 
c  g  C  ci  C  and  the  type  of  search  platfonn  p  e  TYPE  c:  P  (i.e.  MPA,  Seaweb,  frigate 

etc).  We  call  the  cardinally  of  the  set  of  cells  searched  the  size  of  the  mission.  Appendix 
A  lists  the  types  of  search  platforms,  their  detection  capabilities,  classification 
capabilities  and  the  maximum  size  of  mission  each  type  can  perform.  The  geographic 
models  in  Appendix  B  describe  how  detection  and  classification  capabilities  vary  by  cell. 
The  defender  may  have  several  platforms  of  the  same  type  (e.g.  two  frigates  and  three 
Marine  Patrol  Aircraft)  but  each  platform  is  only  capable  of  performing  the  missions 
m  g  M p  c  M  generated  for  that  type  of  platfonn. 

We  begin  generating  missions  by  building  all  the  possible  shapes  for  the  set  of 
cells  in  the  mission  for  each  mission  size.  We  do  not  attach  cells  in  a  mission  diagonally 
since  the  attacker  can  pass  through  such  an  attachment  without  penalty.  We  allow 
rectangle  shapes  of  all  sizes,  “L”  shapes  of  size  three,  and  stair-step  shapes  of  even  sizes 
(Figures  8-10). 
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Figure  8.  Generated  mission  shapes  of  size  three. 


10 


2 

4 

3 

5 

6 

1 

Figure  9.  Generated  Mission  Shapes  of  Size  Six. 
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Figure  10.  Some  Restricted  Mission  Shapes  of  Size  Six.  Arrows  Depict  Zero  Penalty 

Paths. 

Once  we  have  generated  all  the  shapes  of  each  size,  we  generate  all  the  sets  of 
search  platform  navigable  cells  that  can  be  searched  by  each  shape.  These  sets  of  cells 
cannot  include  non-navigable  cells  or  entry  cells.  We  then  combine  these  sets  of  cells 
with  our  platform  types  to  create  our  candidate  missions.  As  discussed  in  Appendix  A, 
some  searcher  types  are  able  to  perform  missions  of  larger  size  than  others.  For  example, 
a  frigate  is  capable  of  searching  missions  up  to  size  three,  while  a  MPA  can  search  up  to 
size  ten. 

We  calculate  the  log  evasion  probability  for  every  cell  in  a  mission  using  both 
detection  probability  and  the  classification  probability  classPiC.  The  data  and  assumptions 
used  to  determine  classification  probability  are  detailed  in  Appendices  A  and  B.  We  use 
the  random  search  formula  to  calculate  the  probability  of  detection.  This  results  in  a 
conservative  estimate  of  the  probability  of  detection.  With  w  being  the  sensor  sweep 
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width,  v  being  speed  and  t  being  the  time  the  attacker  can  transit  a  cell,  and  A  being  the 
area  cell  being  searched,  then  by  the  random  search  formula  (Wagner  1999  p.174): 

Pd  =  1  _e~wvt!  A  ' 

For  moving  search  platforms  such  as  MPA  we  use  the  searcher  speed  for  v.  For 
fixed  platforms  such  as  Seaweb  we  use  the  attacker’s  speed.  We  call  wvt  /  A  the 
coverage  factor  cfp,m,c ■  The  calculation  of  coverage  factor  for  every  platform  type  is 
detailed  in  Appendix  A.  Additional  data  used  in  the  calculation  of  coverage  factor  is 
contained  in  Appendix  B.  Because  the  log  evasion  probability  is  lower  when  the  attacker 
transits  diagonally  than  rectilinearly  (horizontally  or  vertically),  we  calculate  separate 
diagonal  and  rectilinear  and  log  evasion  probabilities. 


Indices  and  index  sets  [-cardinality) 

ceC  cells  with  horizontal,  vertical  coordinates  (alias  cl,  c2)  [—500] 

meM  possible  missions  [-50,000] 

p  eP  defending  platforms  [-10] 

p  £  TYPE  cc  P  type  of  defending  platform  [5] 

m  e  M p  ci  M  missions  platform p  can  perform  [-1 5,000] 

c  £  C  c:  C  cells  patrolled  by  platform  p  while  performing  mission  m  [-20] 


Data  [units] 

cfpim,c  coverage  factor  of  platform  p,  performing  mission  m,  in  cell  c  (See 

Appendices  A  and  B)  [nondimensional]. 

classpc  classification  probability  given  detection  of  platform  p  in  cell  c  (See 

Appendices  A  and  B)  [probability] 


size,,,  number  of  cells  in  mission  m  [1-20] 
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Calculated  Data  [units] 


ev  _  recp  m  c  log  likelihood  that  an  attacker  traversing  cell  c  rectilinearly  would  evade 

detection  and  classification  by  defender  p  performing  mission  in  patrolling 
cell  c  [log  likelihood] 

ev_recpmc  =  log[(l  -(1  -e  Cfp’m’c)classpc\  Vp  eP,  Vce  Cpm,Mm  e  Mp 

evdiag  log  likelihood  that  an  attacker  traversing  cell  c  diagonally  would  evade 

detection  and  classification  by  defender  p  performing  mission  m  patrolling 
cell  c  [log  likelihood] 

ev_diagpmc  =  log[(l - (1  - e^cfp'm’c)classpc]  Vp  ef,  Vce  Cpm,Wm  e  Mp 

Once  we  have  these  log  evasion  probabilities,  we  can  calculate  the  log  probability 
that  the  attacker  traversing  a  cell  will  evade  detection  and  classification  given  an 
incumbent  defender  solution.  We  need  to  account  for  cases  where  the  attacker  traverses  a 
cell  rectilinearly  and  diagonally. 

New  Data  [units] 

% p,m  Incumbent  solution.  1  if  platfonn  p  performs  mission  m,  0  otherwise 

[binary] 

New  Calculated  Data 

e\r_recX c  log  of  probability  that  an  attacker  will  evade  detection  and  classification 

traversing  cell  c  rectilinearly  [log  likelihood] 

recX  c  =  Yj  ev  _  recpmc  S'  p,m  VceC 

peP,meM  p 

^v^_diagX c  log  of  probability  that  an  attacker  will  evade  detection  classification 
traversing  cell  c  rectilinearly  [log  likelihood] 

diagX  c  =  ^  ev  _  diagpmc  X  P,m  VceC 

peP,meM  p 
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c. 


CELL  DATA  TO  ARC  DATA  CONVERSION 


We  now  have  cell-based  log  evasion  probabilities  for  every  cell  in  the  network. 
We  need  to  covert  these  to  arc-based  log  evasion  probabilities.  Figure  11  depicts  an 
example  path  from  cell  (i2,  j  1)  to  cell  (il,  j3).  As  the  attacker  uses  arc  [(i2,  j  1),  (il,  j2)], 
he  is  subject  to  the  diagonal  log  evasion  of  cell  (i2,  jl)  for  half  the  arc  and  the  diagonal 
log  evasion  of  cell  (il,  j2)  for  half  the  arc.  Similarly,  as  the  attacker  uses  arc  [(il,  j2),  (il, 
j3)]  he  is  subject  to  the  rectilinear  log  evasion  of  cell  (i2,  jl)  for  half  the  arc  and  the 
rectilinear  log  evasion  of  cell  (il,  j3)  for  half  the  arc.  If  we  assume  independence,  we  can 
add  x/i  the  log  evasion  of  the  tail  cell  to  Vi  the  log  evasion  of  the  head  cell  to  calculate  the 
total  log  evasion  for  the  arc. 


jl 

j2 

j3 

il 

V 

■ 

Log  evasion  cell  (i2,  j  1) 

i2 

• 

\  ; 

--* 

Log  evasion  cell  (il ,  j2) 

Log  evasion  cell  (il ,  j3) 

Figure  11.  Example  half  arc  log  evasion  calculation.  The  total  log  evasion  for  the  path 
equals  x/2  the  diagonal  log  evasion  for  cell  (i2,  j  1)  plus  A  the  diagonal  log  evasion 
for  cell  (il,  j2)  plus  V2  the  rectilinear  log  evasion  for  cell  (il,  j2)  plus  Vi  the 
rectilinear  log  evasion  for  cell  (il,  j3). 


New  indices  and  index  sets  f—cardi nal itvl 
d  e  Dcl  c2  =  D  cell  adjacencies,  or  traversal  arcs  [-4,000] 

d  e  Drec  rectilinear  arcs 

d  e  Ddi  diagonal  arcs 

New  Calculated  Data  [units] 

evX d  log  of  probability  that  an  attacker  will  evade  detection  and  classification 


traversing  arc  d  [log  likelihood] 
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&X,t=\<y_recXrl+&_recXc2)  WeDm 

^Xd  =  \<^- dia§xci  +  diagX c2)  \Jd  e  Ddiag 


ev 


d,p,m 


log  likelihood  that  an  attacker  traversing  arc  d  would  evade  detection  and 


classification  by  defender  p  performing  mission  m  [log  likelihood] 


evd.n,m  =  T (ev  —  ™cD.m.cl  +ev- recD.m.ci )  V d  e  Drec ,meM 


evd,P,m  =^ev_diagpmcX+ev_diagpmc2)  Vd  e Ddiagjn  e Mp 


The  log  evasion  on  an  arc  is  more  negative  when  the  probability  of  detection  and 
classification  is  high.  On  arcs  where  there  are  no  defenders  present  the  log  evasion  is  at 
its  maximum  value  of  zero  (i.e.,  evasion  probability  one).  Thus,  the  attacker  incurs  no 
penalty  for  traversing  arcs  without  a  defender  present.  This  allows  the  attacker  in  an 
evasion-maximizing  model  to  take  arbitrary  and  unrealistic  meandering  paths,  cycles,  and 
sub-tours.  The  defender  will  attempt  to  interdict  these  meandering  paths  with  an 
unrealistic  sensor  platform  laydown  rather  than  the  more  reasonable  and  realistic  paths 
between  an  entry  cell  and  a  goal  cell.  To  prevent  the  attacker  from  taking  meandering 
paths  we  add  a  small  penalty  to  every  arc  called  fuel.  If  fuel  is  kept  large,  it  induces  no 
unreasonable  behavior  on  the  part  of  the  attacker.  In  probability,  this  represents  a  small 
background  chance  that  an  attacker  will  be  detected  by  complete  accident,  even  though 
the  detection  is  not  made  in  a  cell  with  an  assigned  search  mission.  In  the  absence  of  a 
defender  the  attacker  will  take  the  shortest  path  to  a  goal  cell,  and  the  attacker  will  take  a 
longer  path  to  a  goal  cell  in  order  to  avoid  a  defender. 

New  Data  [units] 

fuel  fuel  penalty  for  traversing  a  single  arc  [log  likelihood  of  evasion] 

Updated  Calculation 

9vX d  =  fuel  +  ^  (ft?_  recXcl  +  recX c2 )  V d  e  Drec 
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EW„  =  fueld-i&diagX,  ^_diagX,2)  W  E  Dt,v 

D.  MODEL  FORMULATION 

Once  we  have  both  the  cell-based  and  arc -based  log  evasion  data,  we  use  a 
modification  of  the  formulation  used  by  Ghaffar  in  his  planning  to  optimize  sensor 
placement  against  terrorist  “go-fast”  boats  (Abdul-Ghaffar,  2008,  pp.  14-19). 

1,  The  Attacker 

The  attacker  has  a  single  SPSS  that  can  choose  to  enter  a  network  at  any  one  of  a 
number  of  entry  cells  c  e  E ,  traverse  a  set  of  cell-to-cell  arcs  d  e  D  to  reach  (and  exit 
the  network  at)  any  one  of  a  number  of  goal  cells  c  eG  .  Traversing  each  arc  carries  a 
risk  of  detection  and  classification  the  attacker  cannot  control,  and  the  log  likelihood  that 

an  attacker  will  evade  detection  and  classification  while  traversing  arc  d  is  evX d .  The 
attacker  seeks  paths  that  maximize  the  log  likelihood  of  evading  detection. 

We  express  the  attackers’  planning  problem  with  the  model  SSV_SP(evV)  : 

New  Indexes  and  index  sets  [-cardinality] 

cefcC  cells  where  an  attacker  can  enter  the  network  [~10] 
ceGcC  goal  cells  [—10] 

Variables  [unitsl 

ENTERc  1  if  the  attacker  enters  network  at  entry  cell  c,  0  otherwise  [binary] 

Yd  1  if  the  attacker  traverses  arc  d,  0  otherwise  [binary] 

GOAL.  1  if  the  attacker  exits  network  at  goal  cell  c,  0  otherwise  [binary] 

Formulation  [dual  variables] 
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zmax (^X)  =  max  X  dYd  (AO) 

deD 

s.t.  YjENTERc  -  1  [a]  (Al) 

cgE 

I  I  r, 

dG°c,c2  deDcl,c 

-ENTERc\cgE+GOALc\cgG<0  VceC  [/?]  (A2) 

-^G044<-1  [£]  (A3) 

ceG 

Discussion 

The  attackers’  objective  (AO)  is  to  maximize  the  total  expected  log  likelihood  that 
attackers  traversing  arcs  from  entry  cells  on  paths  to  goal  cells  evade  detection  and 
classification  (or,  equivalently,  to  maximize  the  joint  probability  that  they  evade 
detection  and  classification  over  all  the  arcs  they  choose  to  traverse).  Constraint  (Al) 
limits  the  entry  into  the  network  to  the  entry  cells,  each  constraint  (A2)  forces 
conservation  of  flow  at  a  cell  in  the  network,  and  constraint  (A3)  limits  the  exit  from  the 
network  to  goal  cells.  This  linear  program  will  produce  an  intrinsically  integral  solution 
Y*  (e.g.,  Ahuja  et  al.  1993,  pp.  447-449).  (The  Greek  notation  associated  with  each 
constraint  denotes  a  linear  programming  dual  variable  used  later.) 

2.  The  Defender 

The  defender  controls  a  set  of  search  platforms  (e.g.,  MPA,  frigates,  etc.)  p  e  P 
that  may  each  be  performing  a  mission  m  e  M  .  The  log  likelihood  that  an  attacker 

traversing  arc  d  evades  detection  and  classification  by  defender  platform  p  performing 
mission  m  is  evd  m  .  The  defender  seeks  positions  for  his  search  platforms  to  minimize 

the  total  log  likelihood  of  attackers  evading  his  surveillance.  We  express  the  defender’s 
problem  as  follows  DMIN(Y) 

New  indices  and  index  sets  f—cardi nal ilyl 

m  e  M  possible  missions  [-50,000] 

p  <eP  defending  platforms  [-10] 
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m  e  M  cz  M  missions  platfonn  p  can  perfonn  [-10,000] 


New  data  [units] 


ev 


d,p,m 


Variables  [units] 


Z 


log  likelihood  that  an  attacker  traversing  arc  d  would  evade 

detection  and  classification  by  defender  p  perfonning  mission  m 
[log  likelihood] 

1  if  attacker  traverses  arc  d,  0  otherwise  [binary] 

=  1  if  platform  p  performing  mission  in  in  cell  c,  0  otherwise 
[binary] 

total  log  likelihood  of  evading  detection  and  classification 
[log  likelihood] 


Formulation 


Zmin 

(Y)  =  min  Z 

^,z 

(DO) 

s.t. 

Z  >  V  ev. 

d,p,m 

deD,peP 

meMp 

Y  Y 

'd^p.m 

(Dl) 

s.t. 

Z  xp,.  s  1 

meMp 

\/p  e  P 

(D2) 

Z  xp,.  5 1 

psP 

Vm  e  M p 

(D3) 

Z  XPP 

( p p  acgC pm 

,  <  1  VceC 

(D4) 

^{0,1} 

Vp  e  P,  Vm  e  M p 

(D5) 

Discussion 

(DO)  introduces  the  objective,  and  constraint  (Dl)  defines  the  objective  variable 
as  the  minimum  upper  bound  on  total  log  likelihood  of  evasion.  Each  constraint  (D2) 
allows  a  defender  platform  to  perfonn  at  most  one  mission,  each  constraint  (D3)  allows  a 
mission  to  be  perfonned  by  at  most  one  platform,  each  constraint  (D4)  allows  at  most  one 
platform  occupy  a  cell,  and  (D5)  stipulates  a  binary  decision  for  each  platfonn  mission. 
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3. 


Defender- Attacker  (D-A)  Model 


We  initially  assume  that  attacker  has  perfect  knowledge  of  these  mission  choices 
and  responds  by  choosing  a  path  that  maximizes  the  cumulative  log  evasion  probability 
created  by  these  mission  choices.  These  sequential  actions  by  the  defender  and  the 
attacker  are  a  type  of  Stackelberg  game  represented  as  a  defender-attacker  (D-A)  model 
(e.g.,  Brown  et  ah,  2006).  The  following  initial  formulation  SPSS-MINMAX  is  a  min- 
max  formulation  that  cannot  be  solved  through  conventional  means. 

Z  =  min  max  V  ev.  n  X  n  Y, 

2  x  Y  d,p,m  p,m  d 

’  deD, 

peP,meM  „ 

s.t.  (Al)  -  (A3)  and  (D 1)  -  (D4) 

4.  Dual  Integer  Linear  Program  Formulation 

A  standard  technique  in  solving  D-A  is  to  temporarily  fix  Z  and  X,  take  the  dual 
of  the  remaining  problem  ( SSV_SP(evX)  ),  and  then  release  Z  and  X  to  obtain  a  single 
integer  linear  program  that  can  be  solved  conventionally  (e.g.,  Brown  et  al.,  2006).  We 
express  this  integer  linear  program  as  SPSS-ILP. 


min 

a -8 

(TO) 

a, (3,8, 
X 

a-  Pc>  0 

VceT 

(Tl) 

S.t. 

(3 c\  ft  cl  ~  Z  evd,p,mxp,m 

e  DcU2 

(T2) 

peP, 

CECpm 

-S  +  P„  >0 

VceG 

(T3) 

Z  Xr:  S  1 

\/p  eP 

(T4) 

meMp 

Z  s  i 

Vc  G  C 

(T5) 

(p,m 

i):niGM  p  acgC pm 

a  >  0 

(3C>  0 
^>0 

^€{0,1} 

VceC 

\/p  e  P,m  e 

Mp{  T6) 
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SPSS-ILP  solves  the  D-A  problem  by  choosing  platform  missions  X,  and  then  recovers 
optimal  attacker  paths  by  fixing  X  and  solving  SSV_SP(&AQ  . 

5.  Decomposition 


The  SPSS-ILP  is  difficult  to  solve  on  a  large  scale.  A  Benders  decomposition 
applies  naturally  to  D-A  problems.  The  standard  Benders  method  takes  the  dual  of 
SPSS-ILP  with  X  fixed  which  leads  back  to  SPSS-MINMAX  (e.g.,  Brown  et  al.  2006). 

We  modify  DMIN(Y)  replacing  equation  (Dl)  with  (D1K).  This  is  the  Bender’s 

decomposition  Master  Problem.  SSV_SP(&AQ  is  the  Bender’s  decomposition 

Subproblem. 

New  index 

k  e  K  decomposition  iteration 

New  Data 

Yk  attacker  plans  for  iteration  k 

DMIN(Y)  Modification 


AY)  =  min  Z 

z,x 


s.t. 


Z> 


Eev  YkX 

e  yd,p,m1d  ^  p,m 

dsD, 

peP,meM  „ 


k=l,...,K  (D1K) 


Constraints  (D2)-(D4)  unchanged. 


The  complete  decomposition  algorithm  is  as  follows: 

Algorithm  SPSS  Decomposition 

Input:  Data  for  defense  problem,  optimality  tolerance  s  >  0  ; 

Output:  {.'-optimal  defender  location  plan  X* ,  and  responding  attacker  plan 

Y*; 
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1 .  Initialize  best  upper  bound  Zm  <—  co  ,  best  lower  bound  ZLB  < - oo , 

define  the  incumbent,  null  defender  plan  X*  <—  X1  <—  0  as  the  best 
found  so  far,  and  set  iteration  counter  K  <—  1 ; 

2.  Subproblem:  Calculate  using  XA.  Solve  subproblem 

SSV_SP(@vX)  to  determine  the  optimal  attack  plan  YA  given  XA ; 
the  bound  on  the  associated  objective  is  Zmax(XA ) ; 

3.  If  (  ZUB  >  Zmax(XA)  )  set  ZUB  <—  Zmax (XK  )  and  record  improved 
incumbent  defender  plan  X*  <—  XA ,  and  responding  attacker  plan 

Y*  <—  Ya  ; 

4.  If  (ZUB  ~ZLB  <  s )  go  to  End; 

5.  Master  Problem:  Given  attack  plans  Y*,  k=l,...K,  attempt  to  solve 
master  problem  DMIN(Y)  to  determine  an  optimal  defender  plan 
XA+1 .  The  bound  on  the  objective  is  Zmin  (Y) ; 

6.  If  ZLB<Zmm( Y)  set  ZLB  <—  Zmin(Y) ; 

7.  If  ( ZUB  -  ZLB  <  s  )  go  to  End; 

8.  Set  K  *—  K +1  and  go  to  step  (2)  (Subproblem); 

9.  End:  X*  is  an  ^-optimal  defender  solution,  and  Y*  is  the  attacker 
response  to  that  plan. 

6.  Binary  Relaxation 

We  initially  assume  that  the  defender’s  platforms  are  operating  overtly  and  the 
attacker  has  perfect  knowledge  of  (“observes”)  defender  actions. 
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The  assumption  of  perfect  attacker  knowledge  is  a  conservative  assumption,  and 
perhaps  not  a  very  realistic  one  in  light  of  the  low-technology  nature  of  an  SPSS.  In  the 
covert  case,  we  assume  that  the  attacker  knows  the  number  and  type  of  platfonns  that 
may  be  searching  but  does  not  know  their  exact  disposition.  In  this  case,  the  optimization 
becomes  a  simultaneous  two-person  zero  sum  game  (Von  Neumann  et  ah,  2004).  The 
attacker  may  need  to  randomize  path(s)  and  the  defender  may  need  to  randomize  mission 
assignments. 

The  formulations  of  the  overt  case  and  covert  case  are  nearly  identical.  In  order 
to  formulate  the  simultaneous-play  game  in  the  covert  case  the  binary  restriction  of  X  and 
Y  are  relaxed  to  continuous  restrictions  (Brown  et  ah,  2008,  p.  97).  The  resulting  X  is  the 
probability  that  the  defender  chooses  to  employ  a  platform  in  a  mission.  We  refer  to  this 
probability  as  a  mixed  strategy,  though  it  may  turn  out  in  practice  to  by  a  pure  one  (i.e., 
all  probability  devoted  to  just  one  mission). 

E.  SOFTWARE  IMPLEMENTATION 

We  implement  the  SPSS  decomposition  algorithm  using  General  Algebraic 
Modeling  System  (GAMS)  (GAMS,  2009)  with  ILOG  CPLEX  (ILOG,  2007)  set  as  the 
linear  programming  solver.  We  use  Visual  Basic  for  Applications  (VBA)  to  support  this 
GAMS  code.  We  calculate  ev_rec  and  ev _diagpmc  in  VBA  using  the  assumptions 

from  Appendix  A  with  the  geographic  models  from  Appendix  B  contained  in  Microsoft 
Excel  worksheets.  Additionally,  VBA  uses  the  GAMS  output  to  generate  the  probability 
field  and  mixed  strategy  figures  discussed  in  the  next  chapter. 
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III.  RESULTS  AND  ANALYSIS 


Typical  solutions  of  the  SPSS  decomposition  algorithm  take  5  hours  to  complete 
270  iterations  using  a  windows  based  personal  computer  with  a  Xeon  X5460  processor. 

A  typical  final  iteration  instance  of  SSV_SP(&AQ  has  3,000  variables  and  450 
constraints  with  10,000  non-zero  elements.  A  typical  final  iteration  instance  of 
DMIN(Y)  has  280,000  variables  and  700  constraints  with  4.1  million  non-zero  elements. 

For  our  scenarios,  we  use  five  different  representative  search  platfonn  types.  The 
characteristics  of  these  platfonn  types  are  summarized  in  Table  1  and  detailed  in 
Appendix  A.  The  coverage  factor  cfp,m,c  used  in  the  calculation  of  the  log  evasion  is 
proportional  to  speed  of  the  platform  (or  speed  of  the  SPSS  in  the  case  of  Seaweb)  and 
sweep  width.  The  sweep  width  of  each  sensor  type  varies  by  geographic  factors 
summarized  in  Table  1  and  detailed  in  Appendix  B.  For  each  platform  type,  the 
classification  probability  c/assPiC  used  the  calculation  of  log  evasion  varies  by  the  density 
of  other  shipping  activity  in  each  cell. 


Platform 

type 

max 

size 

Classification 

Probability 

Speed 

(knots) 

Sweep 

width(nm) 

Sweep  width  depends 

on 

MPA 

10 

0.95-0.80 

180 

0-17.4 

wind  speed 

AEW 

10 

0. 5-0.1 

180 

0-43.5 

wind  speed 

Frigate 

3 

0.95-0.80 

10 

10-30 

wind  speed,  water  depth, 

shipping  density 

Submarine 

3 

0. 5-0.1 

10 

10-30 

wind  speed,  water  depth, 

shipping  density 

Seaweb 

1 

(per 

sensor) 

0.5-0. 1 

6 

(SPSS 

speed) 

0-37 

(per  sensor) 

water  depth 

Table  1.  Sensor  platform  characteristics.  For  instance,  a  frigate  can  cover  missions 

searching  up  to  3  cells,  with  classification  probability  0.95-0.80,  with  speed  10 
knots  and  sweep  width  10-30  nautical  miles,  depending  on  factors  including  wind 
speed,  water  depth,  and  shipping  density  in  the  searched  cells. 
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From  Table  1  we  see  that  AEW  aircraft  are  better  at  detection  than  MPA  but  the 
MPA  is  better  at  classification.  Given  its  wider  sweep  width  the  AEW  aircraft’s 
detection  probability  degrades  less  than  the  MPA’s  detection  probability  as  the  mission 
size  is  increased  to  its  maximum  of  10  cells.  We  therefore  expect  the  defender  to  choose 
missions  of  larger  size  for  AEW  aircraft  than  for  MPA.  MPA’s  classification 
performance  does  degrade  with  shipping  density  by  not  to  the  same  degree  as  AEW.  We 
therefore  expect  the  defender  to  use  MPA  in  areas  of  high  shipping  density  compared  to 
AEW.  Both  of  these  platforms  are  best  used  in  areas  where  the  wind  speed  is  low.  The 
detection  performance  of  both  the  frigate  and  the  submarine  are  hampered  by  their  slow 
speed.  The  detection  performance  of  these  platforms  degrades  quickly  as  the  size  of  the 
mission  increases.  We  expect  most  selected  frigate  and  submarine  missions  to  be  of  size 
one.  Though  both  platforms  are  affected  by  acoustic  conditions  (wind,  water  depth,  and 
shipping  density),  the  submarine’s  classification  probability  is  affected  more  acutely  by 
shipping  density  than  the  frigate.  We  expect  the  defender  to  avoid  using  the  submarine  in 
areas  of  high  shipping  density.  Water  depth  has  the  most  influence  on  Seaweb  detection 
performance.  In  deep  water,  a  ten-sensor  Seaweb  network  has  similar  detection 
performance  to  an  MPA.  Like  the  submarine  and  AEW  aircraft,  Seaweb  is  relatively 
poor  classifier  and  is  better  used  in  areas  of  low  shipping  density. 

We  restrict  the  total  mixed  strategy  probabilities  of  all  platform  assignments  to 
any  given  cell  be  less  than  or  equal  to  one.  Sending  a  platform  to  a  given  cell  with  too 
high  a  probability  precludes  the  use  of  that  cell  for  other  platfonns.  A  relatively  poor 
sensor  platform,  such  as  the  submarine  in  a  desirable  cell,  might  keep  that  cell  from  being 
used  by  more  capable  platfonns.  These  less  capable  platforms  can  be  useful  in  confined 
areas,  such  as  chokepoints.  It  might  be  better  to  employ  a  frigate  in  a  chokepoint  and 
allow  an  MPA  to  spread  its  effort  of  over  more  cells  in  less  confined  areas. 

A.  CARIBBEAN  BASELINE  SCENARIO 

We  choose  the  Caribbean  as  the  setting  for  our  first  scenario  (Figure  12).  The 
geography  of  the  Caribbean  is  more  conducive  to  model  validation  than  the  Pacific. 
There  are  several  chokepoints  that  the  SPSS  must  transit  to  reach  any  goal.  Given  a 
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reasonable  number  of  entry  and  goal  cells,  a  valid  solution  should  make  use  of  these 
chokepoints.  The  Caribbean  is  also  a  plausible  location  for  export  of  SPSS  technology. 
Though  SPSS  have  operated  nearly  exclusively  in  the  Eastern  Pacific,  the  Caribbean  has 
been  used  as  a  transit  zone  for  other  vehicles  used  in  drug  smuggling,  such  as  go-fast 
boats  and  aircraft.  For  this  baseline  scenario,  we  use  one  each  of  MPA,  AEW,  frigate  and 
submarine  platfonn  types.  We  assume  that  we  have  one  of  these  assets  on  station  at  all 
times.  In  addition,  we  employ  two  separate  Seaweb  networks  with  10  sensors  each. 


Figure  12.  Caribbean  Baseline  Scenario  Entry  and  Goal  Cells.  Each  cell  is  60  nautical 
miles  square.  This  depiction  shows  latitudes  08-23N  and  longitudes  64-88W.  On 
the  west  is  the  Mexican  and  Central  American  coast,  on  the  south  is  Columbia, 
where  each  “E”  cell  is  a  possible  entry,  and  to  the  north  is  the  Gulf  coast  of  the 
U.S.,  where  each  “G”  cell  is  a  possible  goal.  Islands  present  obstructions  to 
navigation  for  both  attacker  and  defender. 


1.  Caribbean  Mixed  Strategy 

The  mixed  strategy  we  select  is  the  probability  the  defender  should  adopt  each 
mission  for  any  given  time  epoch.  The  number  in  each  cell  indicates  the  probability  a 
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search  platform  is  sent  to  a  cell.  The  higher  probability  cells  are  darkly  shaded,  and  the 
lower  probability  cells  are  lightly  shaded.  Figures  13-17  depict  individual  platform 
mixed  strategies. 


Figure  13.  Caribbean  MPA  Mixed  Strategy.  Cells  with  dark  outlines  are  entry  and  goal 
cells.  Numeric  values  are  search  platform-to-cell  probabilities  of  assignment. 

The  shading  distinguishes  qualitatively  between  relatively  high  search 
probabilities  and  low  ones.  We  have  chosen  a  probabilistic  barricade,  utilizing 
available  geographic  chokepoints  to  navigation.  Suggested  missions  are  generally 

of  size  three  or  four  cells. 


MPA  is  used  primarily  to  cover  the  chokepoints.  When  compared  with  the 
shipping  density  model  in  Appendix  B,  we  see  that  MPA  are  sent  to  areas  of  high 
shipping  density.  We  also  see  from  the  wind  model  in  Appendix  B  that  MPA  avoid  any 
areas  where  the  wind  is  so  high  that  the  sweep  width  is  zero.  The  defender  tends  to 
choose  MPA  missions  of  size  three  or  four. 
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Figure  14.  Caribbean  AEW  Mixed  Strategy.  Here  we  see  much  wider-ranging  mission 

assignments. 

AEW  aircraft  tend  to  spread  in  areas  of  open  ocean  but  still  use  these  open  ocean 
areas  to  cover  the  chokepoints.  Most  suggested  missions  are  of  size  nine  or  ten. 
Checking  the  shipping  density  and  wind  models  in  Appendix  B  reveals  that  the  defender 
AEW  aircraft  avoid  areas  with  high  shipping  density  and  high  winds  with  a  few 
exceptions  that  merely  allow  assignment  of  attractive  larger  mission  sizes. 
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Figure  15.  Caribbean  Frigate  Mixed  Strategy. 

The  frigate  is  used  in  a  few  constrained  areas,  but  there  is  also  a  tendency  to  avoid 
areas  that  may  be  better  used  by  other  platforms.  When  comparing  to  the  shipping 
density  and  acoustic  condition  models  in  Appendix  B,  we  see  that  all  of  the  cells  chosen 
are  of  moderate  or  high  shipping  density  and  all  have  moderate  to  good  acoustic 
conditions.  All  but  one  of  the  frigate  missions  suggested  for  the  defender  are  of  size  one. 
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Figure  16.  Caribbean  Submarine  Mixed  Strategy. 

When  comparing  to  the  shipping  density  and  acoustic  condition  models  in 
Appendix  B,  we  see  that  all  of  the  cells  chosen  are  of  low  shipping  density  and  all  have 
moderate  to  good  acoustic  conditions.  The  most  likely  missions  chosen  for  the 
submarine  are  of  size  one. 


29 


Figure  17.  Caribbean  Seaweb  Mixed  Strategy. 

Checking  the  depth  model  in  Appendix  B  reveals  advantageous  placement  of 
Seaweb  in  deep  water  to  cover  the  entry  cells  as  well  as  the  eastern  goal  cells.  The 
chokepoints  are  too  shallow  for  Seaweb  to  be  used  effectively.  Examining  the  shipping 
density  model  in  Appendix  B  reveals  that  we  suggest  placing  Seaweb  in  areas  of  high 
shipping  density  only  when  the  water  is  deep  and  the  cell  is  part  of  a  larger  mission  in  a 
desirable  area  for  search.  The  size  of  the  missions  chosen  by  the  defender  vary  by  water 
depth.  The  larger  mission  sizes  are  in  deeper  water. 
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Figure  18.  Overall  Mixed  Strategy.  This  reminds  us  that  the  prior,  platform-by-platform 
mixed  strategies  are,  in  fact,  an  orchestrated,  simultaneous,  coordinated,  joint 
employment  plan  for  all  of  them.  Over  time  this  shows  how  we  suggest  all  the 

platforms  be  assigned. 

Figure  18  depicts  the  overall  mixed  strategy.  This  shows  where  the  defender 
chooses  to  send  all  assets.  A  cell  is  saturated  if  the  total  mixed  strategy  probability  for  all 
defenders  equals  1.0.  We  see  the  defender  saturates  the  easternmost  chokepoint  with 
search  platforms. 

2.  Caribbean  Probability  Field 

We  combine  the  mixed  strategy  with  sensor  performance  to  create  the  probability 
field  in  Figure  19.  More  darkly  shaded  cells  have  a  high  probability  of  detection  and 
classification  that  the  attacker  wishes  to  avoid.  The  numbers  in  the  probability  field  cells 
are  based  on  a  rectilinear  transit.  The  probabilities  based  on  a  diagonal  transit  are  higher 
but  the  relative  strength  of  these  probabilities  from  cell  to  cell  is  similar  for  rectilinear 
and  diagonal  transits. 

The  path  depicted  by  the  arrow  is  the  worst-case  attacker  path  from  the  defender’s 

point  of  view.  The  probability  of  detection  and  classification  given  this  path  is 
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approximately  0.69.  Any  other  path  taken  by  the  attacker  results  in  a  higher  probability 
of  detection  and  classification.  Of  course,  the  attacker  cannot  always  take  this  worst-case 
path.  If  the  attacker  always  used  this  worst  case  path,  the  defender  would  be  able  to 
defend  against  this  one  path  rather  than  all  possible  paths. 


Figure  19.  Caribbean  Probability  Field  and  Worst-Case  Attacker  Path. 


B.  PACIFIC  BASELINE  SCENARIO 

The  Pacific  is  a  more  realistic  setting  than  the  Caribbean  in  that  drug  cartels 
currently  use  large  numbers  SPSSs  here.  There  are  no  chokepoints,  but  there  are  fewer 
entry  cells  than  goals  cells  (Figure  20).  Unless  there  is  a  marked  difference  in 
environmental  conditions  between  the  areas  bordering  entry  and  goal  cells,  we  expect  the 
defender  to  concentrate  his  assets  to  cover  the  area  surrounding  entry  cells.  As  in  the 
Caribbean,  the  defender  has  one  MPA,  one  AEW  aircraft,  one  frigate,  one  submarine,  and 
two  Seaweb  networks  with  ten  sensors  each.  Figures  21-26  depict  the  Pacific  mixed 
strategy. 


32 


Figure  20.  Pacific  Baseline  Scenario  Entry  and  Goal  Cells.  Each  cell  is  a  60  nautical 
mile  square.  At  the  lower  right  (south-east)  is  the  coast  of  Columbia  with  entry 
cells  “E”,  while  at  the  upper  left  (north-west)  is  the  coast  of  Mexico  with  goal 

cells  “G”. 


33 


1. 


Pacific  Mixed  Strategy 


Figure  2 1 .  Pacific  MPA  Mixed  Strategy.  We  see  that  in  this  case  the  suggested 

probabilistic  barrier  covers  the  entry  cells. 


MPA  is  used  primarily  to  cover  the  entry  points.  When  compared  to  the  shipping 
density  model  in  Appendix  B,  we  see  that  MPA  cover  nearly  all  the  high  shipping  density 
cells  in  the  vicinity  of  the  entry  cells.  We  also  see  from  the  wind  model  in  Appendix  B 
that  MPA  avoids  any  areas  where  the  expected  wind  is  20  knots  or  higher.  The  only  high 
shipping  density  cell  in  the  vicinity  of  the  entry  cells  not  covered  by  MPA  has  20-knot 
winds.  Due  to  lower  winds  and  wider  sweep  widths,  we  preferentially  choose  MPA 
missions  of  larger  size  in  the  Pacific. 
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Figure  22.  Pacific  AEW  Mixed  Strategy.  Similar  to  MPA,  AEW  aircraft  concentrate 

their  search  effort  near  the  entry  cells. 


As  with  MPA,  AEW  aircraft  are  better  used  to  cover  the  entry  points.  As  in  the 
Caribbean,  we  prefer  mission  sizes  of  nine  or  ten.  Checking  the  shipping  density  and 
wind  models  in  Appendix  B  reveals  that  the  defender  AEW  aircraft  avoid  areas  with  high 
shipping  density  and  high  winds  with  exception  of  il5,  j25.  Choosing  this  cell  allows 
AEW  defender  to  fit  a  mission  with  larger  mission  size  within  a  relatively  constrained 
navigable  area. 
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Figure  23.  Pacific  Frigate  Mixed  Strategy. 


The  frigate  avoids  areas  that  may  be  better  used  by  other  platforms  and 
preferentially  searches  relatively  constrained  areas.  From  the  acoustic  model  in 
Appendix  B  we  see  that  the  cells  chosen  by  the  defender  have  good  acoustic  conditions. 
Every  frigate  mission  we  suggest  is  of  size  one. 
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Figure  24.  Pacific  Submarine  Mixed  Strategy. 


Checking  the  shipping  density  and  acoustic  condition  models  in  Appendix  B 
reveals  that  all  of  the  cells  suggested  are  of  low  to  moderate  shipping  density  and  have 
moderate  to  good  acoustic  conditions.  The  submarine  missions  chosen  are  of  size  one  or 
two. 
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Figure  25.  Pacific  Seaweb  Mixed  Strategy.  The  overall,  blended  operational  plan 
concentrates  search  effort  near  the  entry  cells.  The  shaded  cells  with  “0” 
probability  have  non-zero  probabilities  too  small  to  display. 


Checking  the  depth  model  in  Appendix  B  reveals  we  select  a  mixed  strategy  for 
Seaweb  that  closely  matches  the  depth  contours.  Unlike  the  other  sensor  platforms,  we 
devote  some  Seaweb  effort  to  cover  areas  in  the  vicinity  of  goal  cells.  These  areas  are  in 
deep  water.  Comparing  with  the  shipping  density  model  in  Appendix  B  shows  we  only 
choose  to  use  Seaweb  in  areas  of  moderate  to  low  shipping  density.  The  size  of  the 
missions  suggested  vary  by  water  depth.  The  larger  mission  sizes  are  in  deep  water. 
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Figure  26.  Pacific  Overall  Mixed  Strategy.  We  suggest  a  mixed  strategy  that  is  a 

blockade  of  entry  cells  in  Columbia. 


We  see  in  the  overall  mixed  strategy  that  we  suggest  a  probabilistic  blockade  in 
the  immediate  vicinity  of  entry  cells. 

2.  Pacific  Probability  Field 

We  combine  the  mixed  strategy  of  Figures  21-26  with  sensor  performance 
calculations  to  create  the  probability  field  in  Figure  27.  The  path  depicted  by  the  arrow  is 
the  worst-case  attacker  path  from  the  defender’s  point  of  view.  The  probability  of 
detection  and  classification  given  this  path  is  approximately  0.86.  This  probability  is 
higher  than  that  we  achieve  in  the  Caribbean  due  to  a  more  favorable  sensor  environment. 
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Figure  27.  Pacific  Probability  Field  and  Worst-Case  Attacker  Path. 


C.  RESOURCES  VERSUS  RISK 

We  now  explore  the  effect  of  removing  MPA  on  the  Baseline  mixed  strategy  in 
the  Caribbean  (Figures  28-32).  With  the  removal  of  MPA  the  defender  adjusts  the  mixed 
strategy  by  using  the  AEW  aircraft  and  Seaweb  to  cover  the  open  ocean  area  between  the 
entry  cells  and  the  two  western  chokepoints.  These  areas  have  low  to  moderate  shipping 
densities  that  allow  AEW  aircraft  and  Seaweb  to  operate  effectively  (see  Appendix  B  for 
shipping  density  and  depth  models).  The  westernmost  chokepoint  with  high  shipping 
densities  now  has  no  direct  coverage.  Similarly  the  defender  deploys  Seaweb  to  cover 
the  eastern  approaches  to  the  goal  cells  by  deploying  in  deep  water  in  columns  j21  and 
j22.  The  frigate  strategy  complements  Seaweb  in  these  columns  by  covering  high 
shipping  density  cells. 
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Figure  28.  AEW  Mixed  Strategy  for  Caribbean  Scenario  with  no  MPA. 


Figure  29.  Frigate  Mixed  Strategy  for  Caribbean  Scenario  with  no  MPA. 
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Figure  30.  Submarine  Mixed  Strategy  for  Caribbean  Scenario  with  no  MPA. 


Figure  3 1 .  Seaweb  Mixed  Strategy  for  Caribbean  Scenario  with  no  MPA. 
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Next,  to  assess  the  contribution  of  each  platfonn  type,  we  successively  remove 
each  platform  type  from  each  region.  The  results  are  summarized  in  Figures  33  and  34. 
In  both  regions,  we  see  that  removing  the  MPA  has  the  biggest  impact,  followed  by 
AEW,  Seaweb,  frigate  and  submarine.  Given  the  MPA’s  effectiveness  at  both  detection 
and  classification  it  is  not  surprising  that  removing  the  MPA  has  the  biggest  impact.  The 
effect  of  removing  the  frigate  has  almost  as  significant  an  impact  as  removing  Seaweb  in 
the  Pacific.  This  gap  in  impact  between  the  frigate  and  Seaweb  is  much  greater  in  the 
Caribbean.  The  availability  and  location  of  deep  water  is  the  likely  reason.  We  see  in 
Appendix  B  that  the  area  adjacent  to  the  entry  cells  in  the  Pacific  is  relatively  shallow. 
This  makes  Seaweb  a  less  effective  detection  sensor  in  an  area  that  would  otherwise  be 
an  attractive  location  for  Seaweb  deployment. 
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Figure  32.  Caribbean  Resources  Versus  Risk.  Each  bar  shows  the  optimal  worst-case 
probability  of  detection  and  classification  for  the  baseline  scenario,  followed  by  a 
sequence  of  restrictions  removing  each  of  the  indicated  platform  types. 


Figure  33.  Pacific  Resources  Versus  Risk.  Each  bar  shows  the  optimal  worst-case 

probability  of  detection  and  classification  for  the  baseline  scenario,  followed  by  a 
sequence  of  restrictions  removing  each  of  the  indicated  platform  types. 
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D 


IMPROVED  INTELLIGENCE 


We  now  modify  the  Pacific  baseline  scenario  to  explore  the  effect  of  improved 
friendly  intelligence  on  attacker  goal  cells.  We  model  this  improved  intelligence  by 
reducing  the  number  of  goal  cells  to  which  the  attacker  can  send  his  SPSS.  As  shown  in 
Figure  35,  the  number  of  goal  cells  is  now  the  same  and  as  the  number  of  entry  cells. 


Figure  34.  Improved  Intelligence  Scenario  Entry  and  Goal  Cells. 
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Figure  35.  Pacific  Improved  Intelligence  Scenario  Overall  Mixed  Strategy.  In  contrast  to 
the  baseline,  having  better  intelligence  that  restricts  the  number  of  goal  cells  leads 
to  suggested  blockades  of  both  entry  and  goal  cells. 


Now  that  we  have  an  equal  number  of  entry  and  goal  cells,  we  expend  significant 
effort  covering  the  areas  surrounding  both  the  entry  and  goal  cells  rather  than 
concentrating  most  of  the  effort  on  the  entry  cells  (Figure  36).  The  improvement  in 
intelligence  improves  probability  of  detection  and  classification  from  0.86  to  0.89.  In  the 
baseline  scenario  the  defender  is  able  to  effectively  cover  the  entry  cells.  Given  this 
effectiveness  the  ability  to  also  cover  the  goal  cells  in  the  improved  intelligence  scenario 
only  yields  a  small  improvement  over  the  baseline  scenario. 
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E  GEOGRAPHIC  RESTRICTION 

We  modify  the  improved  intelligence  scenario  to  restrict  the  use  of  MPA  and 
AEW  aircraft  to  waters  east  of  the  Galapagos  (column  j  14  and  eastward). 


Figure  36.  Pacific  Geographic  Restriction  Scenario  MPA  Mixed  Strategy.  Areas  that 
cannot  be  searched  by  MPA  or  AEW  are  marked  with  an  “X.” 
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Figure  38.  Pacific  Geographic  Restriction  Scenario  Frigate  Mixed  Strategy. 


48 


1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

11 

12 

13 

14 

15 

16 

17 

18 

19 

20 

21 

22 

23 

24 

25 

26 

27 

28 

1 

2 

3 

4 

5 

zr 

6 

0.27 

0.03 

7 

““ 

_ 

8 

9 

_ 

10 

11 

_ 

■ 

12 

0 

_ 

z 

z 

_ 

13 

14 

0.04 

_ 

_ 

_ 

15 

0 

0 

0 

16 

17 

18 

0 

_ 

I 

19 

0.35 

□ 

20 

0.29 

_ 

1 

21 

0 

0.3 

22 

z 

o 

■ 

23 

24 

_ 

_ 

25 

n 

r 

Figure  39.  Pacific  Geographic  Restriction  Scenario  Submarine  Mixed  Strategy. 


Figure  40.  Pacific  Geographic  Restriction  Scenario  Seaweb  Mixed  Strategy. 
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Figure  41.  Pacific  Geographic  Restriction  Scenario  Overall  Mixed  Strategy. 

We  see  from  Figures  3  7 — 42  the  defender  adjusts  to  the  restriction  on  AEW  and 
MPA  employment  by  using  Seaweb  exclusively  in  cells  west  of  the  Galapagos.  The 
defender  also  sends  the  frigate  west  for  most  of  its  mixed  strategy  covering  two  cells  with 
high  shipping  density  (see  Appendix  B).  The  restriction  reduces  probability  of  detection 
and  classification  from  0.89  to  0.88. 
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IV.  CONCLUSIONS  AND  RECOMMENDATIONS 


This  thesis  develops  an  operational  level  mission-planning  aid  to  counter  Self 
Propelled  Semi-Submersibles  (SPSS)  currently  being  used  by  drug  smugglers  in 
Colombia  to  move  large  volumes  of  cocaine  to  the  Pacific  coast  of  Mexico,  and 
potentially  to  the  Caribbean  coast  of  the  United  States.  This  planning  aid  suggests 
optimal  placement  and  disposition  of  friendly  search  platforms  while  considering  the 
intelligent  response  of  the  enemy  SPSS  operators  to  search  efforts.  While  we  agree  that  a 
successful  operation  against  an  SPSS  requires  detection,  classification,  and  interdiction, 
this  planning  aid  addresses  only  detection  and  classification,  arguably  the  more 
technically  challenging  problems. 

The  planning  aid  uses  a  defender-attacker  optimization  model  to  maximize  the 
defender’s  probability  of  successful  detection  and  classification  of  the  SPSS  through  the 
advantageous  disposition  of  these  search  platforms  against  an  intelligent  attacker 
operating  the  semi-submersible.  We  assume  the  attacker  operating  the  SPSS  has 
imperfect  knowledge  of  defender  platfonn  disposition  but  is  aware  that  there  are 
defenders  that  must  be  avoided.  Given  this  assumption,  the  solution  to  the  defender- 
attacker  model  is  a  mixed  (i.e.,  probabilistic)  strategy  for  the  defender  and  a  least-risk 
path  for  the  attacker. 

We  demonstrate  our  defender-attacker  model  with  both  an  Eastern  Pacific  and 
Caribbean  scenarios  using  five  representative  platfonn  types  whose  detection  and 
classification  performance  varies  by  platform,  and  by  geography.  In  all  of  these  cases, 
we  find  that  our  model  prescribes  a  face-valid  defensive  plan;  defenders  take  advantage 
of  geography  by  positioning  at  chokepoints  in  constrained  waterways,  and  they  provide 
coverage  near  attacker  entry  cells  and  goal  cells  in  the  less  geographically-constrained 
scenarios. 

We  evaluate  the  effect  of  removing  each  class  of  defender  platform  on  the 
defender’s  ability  to  detect  and  classify  a  transiting  SPSS.  Given  our  sensor  perfonnance 
assumptions,  the  removal  of  the  Maritime  Patrol  Aircraft  (MPA)  has  a  much  larger  effect 
on  detection  and  classification  than  removing  any  other  platform.  We  also  demonstrate 
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the  value  of  improved  friendly  intelligence:  in  our  example  our  prior  knowledge  allows 
focus  on  fewer  goal  cells.  This  improved  intelligence  yields  a  corresponding 
improvement  on  probability  of  detection  and  classification.  This  also  results  in  the 
defender  re-directing  some  search  resources  away  from  entry  cells  and  toward  the  now 
easier-to-cover  goal  cells.  Finally  we  demonstrate  the  effect  of  a  geographic  restriction 
on  flight  operations  of  the  MPA  and  AEW  aircraft.  The  defender  adjusts  by  sending 
unrestricted  platforms  to  cover  areas  that  the  aircraft  cannot. 

The  model  solves  our  realistic  sized  scenarios  in  a  few  hours  using  commercially 
available  linear  programming  optimization  software.  We  recommend  further  research 
into  developing  a  heuristic  that  would  solve  these  problems  much  faster  without  requiring 
licensed  optimization  software,  and  preferably  on  a  Navy  Marine  Corps  Intranet  (NMCI) 
computer. 

We  attempted  to  run  the  model  assuming  perfect  attacker  knowledge.  The 
resulting  mixed  integer  program  is  very  hard  to  solve.  We  recommend  further  research 
aimed  at  improving  runtime  for  this  assumption.  We  additionally  recommend  research 
into  a  tri-level  model  where  the  attacker  has  perfect  knowledge  of  some  sensor  platform 
operation  plans,  but  imperfect  knowledge  of  others. 

The  model  ignores  interaction  between  platforms.  We  recommend  further 
research  into  a  model  that  improves  combined  performance  when  two  or  more  platforms 
operating  in  close  company. 

This  thesis  contains  all  data  required  for  the  interested  reader  to  reproduce  any  of 
the  experiments  reported. 
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APPENDIX  A.  SENSOR  PLATFORM  PERFORMANCE 
ASSUMPTIONS  AND  DATA 


This  appendix  presents  the  assumptions  and  data  values  used  to  determine  both 
probability  of  detection  and  probability  of  classification.  These  values  will  vary  for  each 
platform  given  the  searcher’s  employment  decisions. 

Indices  and  index  sets  [-cardinality! 

c  e  Cj  cells  with  horizontal,  vertical  coordinates  (alias  cl,  c2)  [-500] 

m  e  M  possible  missions  [-50,000] 

p  eP  defending  platfonns  [-10] 

p  e  TYPE  c:  P  type  of  defending  platform  [5] 

m  e  M p  c  M  missions  platform  p  can  perform  [-10,000] 

c  e  C  c:  C  cells  patrolled  by  platform  p  while  performing  mission  m  [-20] 

Data  [units] 

cfpmc  coverage  factor  of  platform  p,  performing  mission  m,  in  cell  c 

[nondimensional] 

clasSp'C  classification  probability  given  detection  of  platform  p  in  cell  c 

[probability] 

depthc  depth  of  a  cell  [mn] 

size m  number  of  cells  in  mission  m  [1-20] 

Wp  C  sweep  width  of  a  platfonn  p  searching  in  cell  c  [mn] 

vp  speed  of  platform  p  [nm/hr] 

n  number  of  sensors  in  a  Seaweb  network  [10] 

A.  MOVING  SEARCH  PLATFORM  TYPES 

The  majority  of  the  search  platfonn  types  p  e  TYPE  c  P  move  within  the  cells 
they  are  assigned  in  a  mission  c  e  C  „cC.  These  platforms  are  considered  to  be 


53 


searching  randomly  within  these  cells.  This  is  a  conservative  assumption.  With  w  being 
the  sensor  sweep  width,  v  being  the  searcher’s  speed  and  t  being  the  time  the  attacker  can 
transit  a  cell,  and  A  being  the  area  cell  being  searched,  then  by  the  random  search  formula 
(Wagner  1999  p.  174): 

Pd=  \-e~wvtlA. 

We  call  wvt  /  A  the  coverage  factor  (c/).  We  assume  that  random  search  sensors  spread 
their  effort  evenly  among  the  cells  they  are  assigned.  Because  all  cells  total  an  area  3600 
nm2,  A  becomes  3600 *sizem  where  sizem=\Cpm  \  is  the  number  of  cells  assigned  in  a 

mission.  The  SPSS  can  travel  through  a  cell  diagonally  or  rectilinearly  (horizontally  or 
vertically).  The  transit  length  is  60  nm  for  an  SPSS  transiting  rectilinearly.  We  assume 
the  SPSS  speed  is  a  constant  6  knots,  so  t  is  10  hours  for  an  SPSS  traveling  rectilinearly. 
Substituting  the  sweep  of  the  platform  in  a  given  cell  wp  c  and  the  speed  of  the  platfonn  vp 
yields  the  following  coverage  factor: 

cfP,m,c  =  wp,cvP  * 10  /  (3600  * sizem) . 

— c/* 

Which  gives  the  probability  of  detection  of  1  -  e  pjn’c  for  a  rectilinear  SPSS  transit.  For 
a  diagonal  transit,  we  simply  multiply  cfp,m,c  by  .  We  assume  random  search 
platforms  ignore  any  detections  outside  their  assigned  search  cells.  Thus,  the  probability 
of  detection  outside  the  cells  being  searched  is  zero. 

1.  Maritime  Patrol  Aircraft  (MPA)  Platform  Type 


We  assume  MPA  to  search  at  v maritime  patrol  aircraft  180  nm/hr.  We  assume  radar  is 
the  primary  search  sensor  for  MPA.  There  are  several  models  of  MPA  employed  by  Joint 
Interagency  Task  Force  (JIATF)  South  with  different  radars  of  varying  performance. 
The  APS-137  is  employed  by  both  USCG  C-130’s  and  some  USN  P-3’s  and  is  chosen  as 
the  representative  MPA  radar.  Table  2  is  used  to  determine  APS-137  sweep  width 
(w. maritime  patrol  aircraft,  c)  (USCG,  2004,  pg.  FI-55). 
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Table  H-27  Sweep  Widths  for  Forward-Looking  Airborne  Radar  (AN/APS- 137) 


16  Nautical  Mile  Radar  Range  Scale 
(Sweep  Width  in  Nautical  Miles) 

Object  Type 

On-scene  Surface  Winds  (krs) 

<5 

to  10 

to  15 

to  20 

to  25 

to  35 

to  45 

to  55 

to  65 

>65 

4  to  10  person  life  raft 

12.1 

8.6 

3.1 

0 

0 

0 

0 

0 

0 

0 

17  to  25  foot  recreational  boat 

13.6 

11.9 

8.2 

2.8 

0 

0 

0 

0 

0 

0 

26  to  35  foot  recreational  boat 

16.6 

16.3 

15.4 

9.5 

3.9 

0 

0 

0 

36  to  50  foot  recreational  boat 

21.0 

20.7 

19.9 

18.9 

17.5 

14.7 

9  8 

3.5 

0 

0 

32  Nautical  Mile  Radar  Range  Scale 
(Sweep  Width  in  Nautical  Miles ) 

Object  Type 

On-scene  Surface  Winds  (kts) 

<  5 

to  10 

to  15 

to  20 

to  25 

to  35 

to  45 

to  55 

to  65 

>  65 

17  to  25  foot  recreational  boat 

17.4 

15.7 

12.0 

66 

0 

0 

0 

0 

0 

0 

26  to  35  foot  recreational  boat 

21.7 

14.9 

9.3 

2.1 

0 

0 

36  to  50  foot  recreational  boat 

28.7 

22.7 

11.5 

3.8 

0 

Table  2.  APS-137  Sweep  Widths  (From  USCG,  2004,  pg.  H-55). 


We  assume  MPA  operating  in  32nm  radar  range  scale  and  the  SPSS  to  have  roughly  the 
same  radar  cross  section  as  a  17-25  foot  recreational  boat.  Surface  winds  will  vary  from 
cell  to  cell  based  upon  the  surface  wind  data  presented  in  Appendix  B. 

We  assume  each  MPA  is  capable  of  employing  visual  and  infra-red  (IR)  sensors 
that  are  organic  to  the  aircraft.  Thus,  we  assign  MPA  a  high  probability  of  classification 
given  detection  without  the  aid  of  other  platforms.  We  assume  each  MPA  has  a 
probability  of  classification  given  detection  (class maritime  patrol  aircraft, c)  of  0.95  in  cells  with 
sparse  shipping  density,  0.9  in  cells  with  moderate  shipping  density  and  0.8  in  cells  with 
high  shipping  density. 

We  assume  MPA  is  able  to  search  up  to  10  cells. 

2.  Airborne  Early  Warning  Aircraft  (AEW)  Aircraft  Platform  Type 

We  assume  AEW  aircraft  to  search  at  vatrbome  early  warning  aircraft  =180  nm/hr.  Radar 
is  the  only  search  sensor  for  AEW.  There  are  several  types  of  AEW  aircraft  employed  by 
Joint  Interagency  Task  Force  (JIATF)  South.  The  APS- 145  is  employed  by  both  U.S. 
Navy  E-2’s  and  Customs  and  Border  Patrol  P-3’s  and  is  chosen  as  the  representative 
AEW  radar.  We  derive  the  APS- 145  sweep  width  ( W airborne  early  warning  aircraft, c)  from  the 
APS-137  sweep  width.  Radar  range  depends  on  a  number  of  factors  including  target 
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radar  cross  section.  APS-145  and  APS-137  have  different  radio  frequencies  and  the  radar 
cross  section  of  a  target  does  vary  by  radio  frequency.  However,  we  assume  that  this 
difference  in  radar  cross  section  is  negligible.  Transmitted  power  of  APS-137  is 
approximately  500  kW  (Friedman,  2006,  pg  210)  and  APS-145  is  approximately  1  mW 
(Forecast  International,  2006).  Based  on  the  relative  size  of  the  radomes  of  MPA  and 
AEW  versions  of  the  P-3  (Jane’s  Aircraft  Upgrades,  2009)  the  APS-145  antenna  area  is 
roughly  20  times  APS-137  antenna  area.  Because  radar  range  varies  by  the  fourth  root  of 
antenna  area  and  transmission  power,  a  factor  of  2.5  is  applied  to  the  APS-137  sweep 
width  to  obtain  the  estimated  APS-145  sweep  width  (Wagner,  1999,  pg  1 13). 

AEW  aircraft  can  only  rely  on  radar  for  classification  and  have  a  relatively  poor 
ability  to  classify  targets.  We  assume  AEW  aircraft  have  a  probability  of  classification 
given  detection  (cl ass  airborne  early  warning  aircraft, c)  of  0.5  in  cells  with  sparse  shipping  density, 
0.3  in  cells  with  moderate  shipping  density  and  0.1  in  cells  with  high  shipping  density. 

We  assume  AEW  aircraft  are  able  to  search  up  to  10  cells. 

3.  Frigate  Platform  Type 

We  assume  frigates  search  at  v frigate  =10  nm/hr.  Passive  towed  array  sonar  is 
assumed  to  be  the  frigate’s  primary  search  sensor.  We  assume  sweep  width  (wfrigate,c)  of 
the  towed  array  is  10  nm  in  cells  with  poor  acoustic  conditions,  15  nm  with  moderate 
acoustic  conditions  and  30  nm  with  good  acoustic  conditions. 

Frigates  are  equipped  with  embarked  helicopters  to  aid  in  classification.  Thus, 
we  assume  the  probability  of  classification  given  detection  (class frigate, c)  is  0.95  in  cells 
with  sparse  shipping  density,  0.9  in  cells  with  moderate  shipping  density,  and  0.8  in  cells 
with  high  shipping  density. 

We  assume  frigates  are  able  to  search  up  to  three  cells. 

4.  Submarine  Platform  Type. 

Submarines  are  assumed  to  be  searching  at  v submarine  =10  nm/hr.  We  assume  the 
hull  mounted  passive  sonar  is  the  submarine’s  primary  search  sensor.  We  assume  sweep 
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width  (w submarine, c)  of  the  towed  array  is  10  nm  in  cells  with  poor  acoustic  conditions,  15 
mu  with  moderate  acoustic  conditions  and  30nm  with  good  acoustic  conditions. 

We  assume  submarines  are  able  to  search  up  to  three  cells. 

Submarines  can  only  rely  on  acoustic  data  for  classification  and  have  a  relatively 
poor  ability  to  classify  targets.  We  assume  submarines  have  a  probability  of 
classification  given  detection  (class submarines)  of  0.5  in  cells  with  sparse  shipping  density, 
0.3  in  cells  with  moderate  shipping  density  and  0.1  in  cells  with  high  shipping  density. 

B.  SEAWEB  PLATFORM  TYPE 

Seaweb  is  a  proposed  autonomous  network  of  acoustic  sensors  that  communicate 
with  each  other  through  acoustic  modems.  Data  from  these  sensors  is  routed  to  a 
gateway  buoy  that  communicates  to  operators  through  radio  or  satellite  signal.  A  Seaweb 
sensor  deployed  in  deep  water  is  expected  to  have  detection  ranges  of  approximately  1 0 
times  the  water  depth  ( 10*depthc )  against  noisy  targets  such  as  SPSS  (Rice,  2008).  We 
use  the  random  search  fonnula  to  calculate  the  probability  of  detection.  This  results  in  a 
conservative  estimate  of  the  probability  of  detection.  Unlike  the  moving  search  sensors 
the  Seaweb  network  is  fixed  within  a  cell  once  it  is  installed.  As  a  result  the  speed  of  the 
SPSS  is  used  (v=6  knots)  in  the  calculation.  As  in  the  moving  search  platforms,  10  hours 
is  used  for  t  for  rectilinear  SPSS  transits.  We  define  n  as  the  total  number  of  sensors  in  a 
Seaweb  network.  We  assume  a  Seaweb  search  sensor’s  sweep  width  is  10*depthc. 
Assuming  the  n  sensors  in  a  Seaweb  network  are  spread  evenly  in  the  size,,,  assigned  cells 
results  in  the  following: 

cf seaweb, m,c  =  n(l0*depthc)*6*  10  /  (3600  *  sizej  . 

This  gives  of  probability  of  detection  of  l-e  ^seaweh’m’c  for  a  rectilinear  SPSS  transit. 
We  assume  Seaweb  sensors  ignore  any  detections  outside  their  assigned  search  cells. 
Thus,  the  probability  of  detection  outside  the  cells  where  Seaweb  sensors  are  deployed  is 
zero.  We  assume  Seaweb  has  a  probability  of  classification  given  detection  (class seawebt() 


57 


of  0.5  in  cells  with  sparse  shipping  density,  0.3  in  cells  with  moderate  shipping  density 
and  0. 1  in  cells  with  high  shipping  density. 

We  assume  the  sensors  are  deployed  in  a  manner  that  allows  communication  with 
a  gateway  buoy.  A  Seaweb  network  can  be  deployed  in  as  many  as  n  cells.  n=  10  for  all 
scenarios  in  this  thesis. 
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APPENDIX  B.  GEOGRAPHIC  MODEL  ASSUMPTIONS  AND  DATA 


This  appendix  presents  the  assumptions  and  data  used  to  determine  how 
probability  of  detection  and  probability  of  classification  vary  by  geographic  location.  We 
divide  seaspace  into  60nm  by  60nm  square  cells.  The  attacker  chooses  a  directed  path 
from  an  advantageous  entry  cell  node  to  an  advantageous  goal  cell  node.  Those  cells 
with  an  “X”  in  Figures  43  and  44  depict  areas  that  are  non-navigable  and  cannot  be  used 
by  the  defender  and  can  only  be  used  by  the  attacker  if  they  are  entry  cells. 


Figure  42.  Useable  cells  in  the  Caribbean.  Cells  marked  “X”  are  non-navigable  and 
cannot  be  used  by  the  defender  and  can  only  be  used  by  the  attacker  if  they  are 

entry  cells. 
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Figure  43.  Useable  cells  in  the  Pacific.  Cells  marked  “X”  are  non-navigable  and  cannot 
be  used  by  the  defender  and  can  only  be  used  by  the  attacker  if  they  are  entry 

cells. 


A.  SHIPPING  DENSITY 

We  assume  shipping  density  influences  probability  of  classification.  It  is  more 

difficult  to  correctly  identify  an  SPSS  if  is  among  other  shipping  than  if  it  is  transiting  an 

area  alone.  Different  types  of  shipping  make  classification  difficult  for  different  types  of 

sensors.  For  example,  sailboats  may  look  very  similar  to  an  SPSS  on  radar  while  fishing 

boats  may  sound  similar  to  an  SPSS  on  acoustic  sensors.  We  ignore  the  effect  of 

different  types  of  shipping  by  assuming  that  the  shipping  densities  we  use  convert  to  an 

equal  number  of  ships  that  might  confuse  radar  or  acoustic  sensors.  Figures  45  and  46 

show  a  one -month  plot  of  ships  that  report  metrological  data  to  the  National  Oceanic  and 

Atmospheric  Administration  (NOAA)  through  the  Shipboard  Environmental  data 
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Acquisition  System  (SEAS)  that  we  use  to  decide  which  cells  have  high,  moderate,  or 
sparse  shipping  densities  (NOAA,  2008).  We  then  assign  a  3  for  high  shipping  density  2 
for  moderate  and  1  for  sparse  in  the  shipping  density  model. 
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Figure  44.  SEAS  Data  and  Corresponding  Caribbean  Shipping  Density  Model.  Densities 
are  scored  with  discrete  values  from  1  (low)  to  3  (high)  (From  NOAA,  2008). 
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Figure  45.  SEAS  Data  and  Corresponding  Pacific  Shipping  Density  Model.  Density 
scores  vary  from  1  (low)  to  3  (high)  (From  NOAA,  2008). 
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B. 


SURFACE  WINDS 


As  discussed  in  Appendix  A,  surface  winds  influence  radar  sweep  widths.  We 
use  a  10-day  average  of  surface  winds  taken  from  Special  Sensor  Microwave/Imager 
(SSM/I)  satellite  data  as  a  representative  wind  distribution  for  both  the  Caribbean  and 
Pacific  (NOAA,  2009). 


Figure  46.  10-Day  Average  Wind  February  06  2009  (From:  NOAA,  2009).. 


We  use  the  10-day  average  winds  from  Figure  47  to  construct  the  wind  models  shown  in 
Figures  48  and  49. 
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Figure  47.  Caribbean  Wind  Model  (average  wind  in  knots). 
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Figure  48.  Pacific  Wind  Model  (average  wind  in  knots). 

C.  WATER  DEPTH 

As  discussed  in  Appendix  A,  water  depth  influences  Seaweb  sensor  performance. 
Water  depth  data  from  Google  Earth  5.0  is  used  to  build  the  water  depth  models  in 
Figures  50  and  51  (Google  Earth,  2009). 
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Figure  49.  Caribbean  Water  Depth  Model  (average  depth  in  hundreds  of  feet). 
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Figure  50.  Pacific  Water  Depth  Model  (average  depth  in  hundreds  of  feet). 


D.  ACOUSTIC  CONDITIONS 

The  performance  of  passive  acoustic  sensors  such  as  a  towed-array  sonar  or  a 
submarine’s  hull-mounted  sonar  depend  on  a  number  of  factors.  We  only  consider  a  few 
of  these  to  help  build  a  representative  distribution  of  cells  with  good,  moderate  and  poor 
acoustic  conditions  for  both  of  our  geographic  areas  of  interest.  The  factors  considered 
are  shipping  density,  wind,  and  water  depth.  Higher  shipping  density  results  in  higher 
ambient  noise  that  degrades  acoustic  sensor  performance.  Higher  winds  lead  to  higher 
sea  states  that  also  increase  the  ambient  noise.  Shallow  water  has  less  favorable  acoustic 
conditions  than  deeper  water.  In  order  to  combine  all  these  factors,  we  assign  a  1,  2  or  3 
for  each  factor.  We  have  already  expressed  shipping  density  in  this  manner.  For  wind, 
we  assign  1  for  average  winds  less  than  or  equal  to  10  knots,  a  2  for  winds  greater  than  10 
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knots  and  less  than  or  equal  to  20  knots,  and  3  for  winds  greater  than  20  knots.  For  water 
depth,  we  assign  1  for  average  water  depth  greater  than  or  equal  to  5,000  feet,  a  2  for 
water  depth  less  than  5,000  and  greater  than  or  equal  to  1,000  feet,  and  3  for  water  depth 
less  than  1,000  feet.  We  add  these  three  factors  to  obtain  the  acoustic  condition  models 
in  Figures  52  and  53.  We  consider  numbers  in  the  range  3-5  to  represent  good  acoustic 
conditions,  6-7  to  represent  moderate  acoustic  conditions,  and  8-9  to  represent  poor 
acoustic  conditions. 


Figure  5 1 .  Caribbean  Acoustic  Condition  Model.  Cell  entries  in  the  range  3-5  represent 

good,  6-7  moderate,  and  8-9  poor  conditions. 
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Figure  52.  Pacific  Acoustic  Condition  Model.  Cell  entries  in  the  range  3-5  represent 

good,  6-7  moderate,  and  8-9  poor  conditions. 
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